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Abstract 

Plasma toroidal metric singularities in helical devices and tokamaks, giving rise to mag- 
netic surfaces inside the plasma devices are investigated in two cases. In the first we con- 
sider the case of a rotational plasma on an helical device with circular cross-section and 
dissipation. In this case singularities are shown to place a Ricci scalar curvature bound on 
the radius of the surface where the Ricci scalar is the contraction of the constant Rieman- 
nian curvature tensor of magnetic surfaces. An upper bound on the initial magnetic field 
in terms of the Ricci scalar is obtained. This last bound may be useful in the engineer- 
ing construction of plasma devices in laboratories. The normal poloidal drift velocity is also 
computed. In the second case a toroidal metric is used to show that there is a relation be- 
tween singularities and the type of tearing instabilities considered in the tokamak. Besides, in 
this case Ricci collineations and Killing symmetries are computed. The pressure is computed 
by applying these constraints to the pressure equations in tokamaks. PACS numbers: 

02.40.Hw-Riemannian geometries 
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I Introduction 



Geometrical techniques have been used with great success [1] in Einstein general relavity have 
been also used in other important areas of physics, such as plasma structures in tokamaks as 
been clear in the seminal book by Mikhailovskii [2] to investigate the tearing and other sort of 
instabilities in confined plasmas [3], where the Riemann metric tensor plays a dynamical role 
interacting with the magnetic field through the magnetohydrodynamical equations (MHD). 
Recentely Garcia de Andrade [4] has applied Riemann metric to investigate magnetic flux 
tubes in superconducting plasmas. Thiffault and Boozer [5] have also applied the methods of 
Riemann geometry in the context of chaotic flows and fast dynamos. In this paper we use 
the tools of Riemannian geometry, also used in other branches of physics as general relativity 
[6], such as Killing symmetries , Riemann and Ricci [7] collineations , shall be applied here to 
generate magnetic nested surfaces in helical devices and tokamaks through the built of Einstein 
spaces obtained from Tokamak plasma metric. This work is motivated by the fact that the 
magnetic surfaces are severely constrained in tokamaks [3] and Killing symmetries are well 
applied every time we have symmetries in the problem as in solutions of Einstein equations 
of general relativity. Equilibrium of these surfaces or their instabilities are fundamental in 
the constructio of tokamaks and other plasma devices such as stellarators where torsion is 
also present. The magnetic surfaces are more easily obtained when symmetries are present. 
This is our main motivation to apply the special Riemann geometrical techniques of Killing 
symmetries and Ricci collineation to obtain magnetic surfaces formed by Einstein spaces. To 
simplify matters we shall consider two usual approximations from plasma physics [3] which are 
the small toroidality or inverse aspect ratio e = -| << 1 where here R represents the external 
radius of the torus and a is its internal radius, and the Shafranov displacement A' << 1 as 
well. We consider examples of two plasma metrics: The first is the metrical of plasma rotation 
in tokamaks where plasma dissipation is taken into account. This metric was used for the 
first time by Tsypin et al [8] to describe dissipative plasmas in the circular cross-section helical 
devices such as HELIACS or Drakon (stellarator) [9]. The pressure of the tokamak is obtained 
from the tokamak Shafranov shift equation. Constant pressure closed ergodic nested surfaces 
in magnetohydrostatics have also been shown by Schief [10] to be generated by solitons. This 
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is another mathematical technique , distinct from ours, is another use ofmathematical theory 
to generate of nested magnetic surfaces in plasmas. The paper is organised as follows: In 
section 2 we review the Ricci tensor techniques and Ricci colineations which are not usually 
familiar to the plasma physicists. In section 3 we solve the Ricci tensor components from the 
plasma metric by considering that nested surfaces are formed by Einstein spaces, where the 
Ricci tensor is proportional to the metric. In section 4 stablishing a geometrical method for 
the classification of tearing instabilities and solve the Ricci collineation equations to find out 
the Killing vectors for the dissipative rotational plasma metric. Conclusions are presented in 
section 5. 



II Ricci collineations from plasma metrics 

Ricci tensor R ik is contructed from the contraction of Riemann tensor RijM in terms of the 
contravariant components of the metric by the expresion 

Ri k = g ll R %J ki (HI) 

and the Ricci scalar R by an extra contraction in terms of the metric tensor as R = g %3 Rij. 
Let us now compute the Riemann space of constant curvature represented by the Riemann 
tensor components 

RijM = A(g ik gji - gug jk ) (II.2) 

where A is a constant which is called de Sitter cosmological constant. Contraction of expression 
(11.2) in two non-consecutive indices, otherwise the symmetry of the Riemann curvature tensor 
Rijki — —RjiH = Rjiik would make them vanish, yields the Einstein space Ricci relation 

Rik = 2Ag tk (II.3) 
The Ricci collineations equations are given by 

[m^rf + R u d k rf + R kl d iV l = (II.4) 

where rf are the components of the Killing vector fj which defines the symmetries of the 
associated space, and di := -£j are the components of the partial derivative operator. This 
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equation is obtained from the more elegant definition in terms of the Lie derivative C n as 

C v R lk = (II.5) 

In the next section we shall construct the magnetic surface as plasma metric singularities and 
in section IV we solve the Ricci collineation equations in terms of the plasma metric above. 



Ill Rotational Plasmas Metric Singularities in Heli- 
cal Devices 

Let us now consider the application of the Tsypin et al metric of a rotational toroidal dissipative 
plasma tokamak, given by the nonvanishing components 

9 - = idw (IIL6) 

& = 4- < iil7 » 

nB 

923 = Rto^- (IH.8) 

913 = -T^tM-^t) ( m -9) 
2nB 2 B 2 

g 33 = R 2 [l - (^-)\[Kcos6}} (111.10) 

7lB 

where the coordinates (a,9,r]) are respectively the internal radius of the torus, which is con- 
stant on the magnetic surface, and the remaining coordinates are the poloidal and toroidal 
angles. The curvature k of the magnetic axis depends in general of toroidal coordinate, and 
the torsion of the magnetic axis is given by r. In this paper a great simplification would 
be possible by considering that the torsion and curvature would be constant. These metric 
allowed Tsypin et al [8] to consider ion viscosity in the plasma. Further ahead we shall also 
consider the computation of the velocity in terms on the singular metric magnetic surface. 
The magnetic field which is chosen in the form B = B [l — r] t cos6 — r] h cos(m6 — nrj)], where 
rj t is the toroidal amplitude of the magnetic field spectrum and rj h is the helical amplitude, 
while (m,n) are the poloidal and toroidal modes of the helical magnetic field. Here we shall 
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consider that the B is constant and that the metric component g 23 — due to to torsion 
weakness assumption. Constancy of B also yields g± 3 also vanishes , which turns the plasma 
helical metric diagonal. Let us now computing the Riemann tensor component 

Rl212 = | &9i2 _ d 2 g u _ d 2 g 22 
dx l dx 2 dx 2 dx 2 dx l dx l 



which yields for the above metric 



#1212 = -^T- (IIL12) 



Since we are looking for nested magnetic surfaces which possess constant curvature the helical 
metric above must satisfy the relation 



which in our case yields 



Rijki = Mdikgji - 9ij9ki) (III. 13) 



#1212 = H9u922 - 912912) (111.14) 



or 



A 

4ttTb 2 

Since the curvature component #1212 is the same in both expressions (111.12) and (111.16), 
equating them results in 



V (111.16) 



By integration one obtains 



4ttB 



Aa 2 

^ = Ta — cT ° ia ° 2 ( IIL17 ) 



where c\ and c 2 are integration constants. Now as usual in general relativity we assume that 
a singularity could be obtained either by the equations gn = 00 and #22 = 0, as happens as 
in Schwarzschild static black hole metric solution of Einstein vacuum field equation. This set 
of equations can be obtained from the helical metric above by setting 0(a) = 0, which from 
(111.17) yields the following second order algebraic equation in the radius a 

Aa 2 



16nB 



+ Cl a + c 2 = (HI.18) 



which solution is 



8nB 0l , 1 , ri A 



a ° A^ [Cl( " 1±[1 "8^ ])] (IIU9) 



Thus we are left with two solutions , the simpler of which is a + = — c\ and 

ao- = [1 - (111.20) 

Note that if ci > solution a is unphysical because the radius of the helical plasma device 
cannot be negative, however the physical solution a ~ yields a magnetic field bound from the 
Ricci scalar 

he 

B » £ id? < ra - 21 > 

To obtain the Ricci collineations we compute first the Ricci scalar components, since the metric 
is diagonal the only surviving components are 

Rn = 9 22 Ri2i2 = ( IIL22 ) 
R22 = 9 U Rm2 = A(p (IH.23) 
These equations yields the Ricci scalar as 

R = g n Rn + g 22 R 22 = A(l + ttB )4> (111.24) 

since in general in plasma devices the field B is very strong we can consider that nB 2 >> B 
which from expressions (111.21) and (111.24) yields an expression which bounds B in terms of 
the Ricci curvature scalar as 

Bo < ^ (III.25) 

This formula can certainly help in the building of new stellarators and toroidal plasma devices 
in general. 

IV Magnetic surface singularities in tokamaks and 
Ricci collineations 

Let us now start by considering the plasma metric given by Zakharov and Shafranov [11] to 
investigate the evolution of equilibrium of toroidal plasmas. The components g ik (i,k=l,2,3) 
and {a,0,z) as coordinates, of their plasma metric are 

g n = 1 - 2A'cos6 + A' 2 (IV.26) 



922 = a 2 (IV.27) 

#33 = (R - A + acosO f (IV.28) 

g 12 = aA'sinO (IV.29) 

where z = asind and the dash represents derivation with respect to a. In our approximation 
the last term in the expression (II. 1) may be dropped. The magnetic surface equations are 
tori with circular cross-section and equations 

r = R- A(a) + acosO (IV.30) 

z = asinO (IV.31) 

We shall consider now just two independent coordinates (x l = a, x 2 = 6) since nested surfaces 
are bidimensional in the case of plasmas, Let us now compute the Riemann tensor components 
in the linear approximation 

Rl212 = | d 2 ffi2 _ d 2 g u _ d 2 g 22 
dx 1 dx 2 dx 2 dx 2 dx 1 dx 1 

Substitution of the plasma metric above into expression (IV. 32) yields the expression 

i?i2i2 = [3A'cos9 - 2] (IV.33) 

It is easy to show that the components i?i 3 i 3 and i? 2 323 both vanishes within our approx- 
imations. At this point we consider that 9 is so small that sinQ vanishes and cosO = 1 
this simplifies extremely our metric and turns it into a diagonal metric where c? 12 = and 
g bb — (gbb)^ 1 (b = 1,2) and this allows us to compute the components of the Ricci tensor 
from the Riemann component. But before that let us compute the use the condition that the 
nested surface is an Einstein space to compute the Riemann component again 

R 1212 = Aa 2 [l - 2A'cos6] (IV.34) 

Since both expressions for the Riemann component i?m2 must coincide, equating expressions 
(IV.33) and (IV.34) yields an expression for the derivative of the Shafranov shift A as 

A' = -6[1 - ^] (IV.35) 
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Integration of this expression yields the value of the shift in terms of the radius a as 

Aa 2 

A = -6a(l + -jy] (IV.36) 

which satisfies the well-known boundary condition A(0) = 0. From these expressions one may 
also compute A" = 2Aa. An important result in plasma physics is that tearing instabilities 
coming from ion or electron currents possess the shift condition A' < 0. This condition would 
be clearly fulfilled from expression unless the A curvature constant would be negative and in 
modulus < —1. Note that this situation is very similar to the condition of favorable or 
unfavorable curvature for the instabilities in plasmas [3]. The main difference is that here 
we are refereeing to Riemann curvature and not to Frenet curvature of the magnetic lines in 
plasmas. This suggests another method to classify geometrically tearing instabilities. Actually, 
since has been shown [12] recently that the Riemann tensor in plasmas can be expressed in 
terms of the Frenet curvature both methods seems to be equivalent. Now let us compute the 
Ricci components Ru and R 22 from the component -R1212 by tensor contraction with metric 
components g 11 and g 22 , which results in the expressions 



and 



which in turn yields the expressions 



and 



Ru = A[l - 2A'] (IV.37) 



R 22 = -[2 + A'] (IV.38) 



diRu = -2AA" (IV.39) 



d l R 22 = -A" (IV.40) 
From equations for i = 1, k = 2 one obtains 

d lV2 = (IV.41) 

d 2 i] x = (IV.42) 

Substitution of these derivatives of the Ricci tensor components into the Ricci collineations 
equations one obtains the following set of PDE equations 

2 [2 + A']^ 1 + 7/ 1 A" = (IV.43) 



which yields 

if = [1 + A'] (IV.44) 
Due to constraint (III. 24) the only solution for the equation 

-[2 + A']d 2 r) 2 - (1 + A') A" = (IV.45) 

is rf = O.This allows us finally to write down the Killing vector as 

77=[(1 + A'),0] (IV.46) 

the last relation which this vector will have to satisfy g^rf) 2 = 1 for the modulus of the 
Killing vector will allow us to determine A and in turn from expression (111.17) will allow us to 
determine the nested surface radius a in terms of the curvature constant A. This implies that 

A' = 1 (IV.47) 

or A = j^a, which satisfies the well-known boundary condition A(0) = 0. Expression (IV. 29) 
yields that A' > which shows physically that the tearing instability cannot come from ion or 
electron currents [2]. Substitution of this result into (111.17) yields 



a = ^A (IV.48) 

The pressure now is easily computed from the expression given in reference 3 

-2nc 2 p^ = -6aJ 2 A - (1 - |-) J 2 (IV.49) 
R R 



where by using the value of A yields 



D 72 

»=3^A* < IV - 5 °> 



one notes that a singularity in the pressure decreases as the curvature constant increases which 
agrees with the reasoning that [11] curvature tends to stabilize the plasma. 



V Conclusions 



In conclusion, we have investigated a method of classification and identification of tearing 
instability, allowing for example to distinguish between tearing instabilities that comes from 

9 



ions and electron currents or not, based on the Riemann curvature constant submanifolds as 
nested surfaces in Einstein spaces. The Killing symmetries are shown also to be very useful 
in the classification of plasma metrics in the same way they were useful in classifying general 
relativistic solutions of Einstein's gravitational equations in four-dimensional spacetime [6]. 
Since as it is well-known [3] the A' behaves as ^ near marginal stability, we must conclude 
that there is a relation between the stability SW > or instability SW < and the positive or 
negative Riemann curvature of the nested surfaces discussed here. Other interesting examples 
of the utility of the is method is the Ricci collineations investigations of the twisted magnetic 
flux tubes and the Arnolds metric for the fast dynamo [13, 14, 15]. Though the examples 
worked here keep some ressemblance to analog gravity models [16] since our black hole analogy 
does not carry much to the metric since the plasma metric we use is a real plasma metric and 
not a pseudo-Riemannian plasma metric built from the scalar wave equation. 

Acknowledgements 

Thanks are due to CNPq and UERJ for financial supports. 



10 



References 

[1] E. Cartan, Riemannian geometry in an orthonormal Frame, (2001) Princeton University 
Press. 

[2] A. Mikhailovskii, Instabilities in a Confined Plasma, (1998) IOP. 

[3] R. White, The theory of toroidally confined Plasmas, revised second edtion (2006) Imperial 
College Press. 

[4] L. C. Garcia de Andrade, Curvature and Torsion effects on carrying currents twisted solar 
loops, (2006) Phys of Plasmas nov issue. 

[5] J. Thiffault and A. H. Boozer, Chaos 11, (2001) 16. L.C. Garcia de Andrade, Physics of 
Plasmas 13, 022309 (2006). 

[6] H. Stephani et al, Exact solutions of Einstein field equations (2003) Cambridge university 
press. 

[7] G. Ricci, Tensor Ana lysis, Boston. 

[8] V.S.Tsypin et al , Phys Plasmas 9 (1998) 3358. V. Tsypin, A.G. Elfimov, R. M. O. 
Galvao.Phys Plasmas 14 (2007). 

[9] M. Wakatani, Stellarator and Heliotron Devices, Oxford (1998). 

[10] W.K. Schief, J. Plasma Physics 69 (2003)465. 

[11] L.E. Zakharov and V. D. Shafranov, Evolution of Equilibrium Toroidal Plasmas in Plasma 
Physics , MIR physics series, Moscow (1981). 

[12] L.C.Garcia de Andrade, Non-Riemannian magnetic flux tubes, Rev Bras Fisica (2006). 

[13] V. Arnold and B. Khesin, Topological Methods in Hydrodynamics, Applied Mathematics 
Sciences 125 (1991) Springer. 

[14] S. Childress and A. Gilbert, Stretch, Twist and Fold: The Fast Dynamo (1996)(Springer). 



11 



[15] B. Khesin, Topology Bounds Energy, in Geometry and Topology in Fluid Flows, Ed. by 
Renzo Ricca (2001) Kluwer Academic Press. 

[16] L.C. Garcia de Andrade, Phys. Rev.D, (2004) 064004. 



12 



